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Abstract 

In  this  work,  a  novel  style  of  straight  flow  channel  tapered  in  height  or  width  is  proposed  to  improve  the  efficiency  of  fuel  utilization  for  PEM  fuel 
cells.  Fuel  channels  of  various  height  and  width  taper  ratios  are  numerically  analyzed  to  understand  their  effects  on  fuel  transport  characteristics 
and  cell  performance.  Influences  of  the  liquid  water  formation  on  the  transport  phenomena  and  cell  performance  are  included  in  the  numerical 
model.  The  present  results  demonstrated  that,  with  the  tapered  channel  designs,  the  flow  area  contraction  along  the  flow  channel  leads  to  increase 
in  fuel  velocity  and  thus  enhances  the  fuel  transport  through  porous  layers,  fuel  utilization,  and  the  capability  of  the  liquid  water  removal.  The 
results  also  reveal  that  the  cell  performance  can  be  improved  by  either  decreasing  height  taper  ratio  or  increasing  width  taper  ratio.  If  the  power 
loss  due  to  pressure  drop  is  not  considered,  the  performance  of  the  fuel  cell  with  the  tapered  flow  channels  is  consistently  improved  with  height 
taper  ratios  decreased  and  width  taper  ratio  increased.  With  the  pressure  loss  considered,  however,  the  best  performance  can  be  obtained  at  the 
height  taper  ratio  ( Ax )  of  0.5  and  the  width  taper  ratio  (Az)  of  1.8  among  the  taper  ratios  studied  in  the  present  work. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

In  PEM  fuel  cells  with  conventional  flow  channels,  reactant 
gases  diffuse  from  channel  inlet  into  catalyst  layer  through  gas 
diffuser  layer  (GDL),  and  then  take  part  in  the  reaction.  Since 
the  diffusion  process  is  not  helpful  for  removing  heat  and  cath¬ 
ode  water,  various  new  styles  of  flow  channel  designs,  such  as 
serpentine,  Z-type  and  interdigitated  designs,  are  proposed  to 
improve  these  matters. 

A  number  of  investigations  were  carried  out  on  serpentine 
flow  channel  designs.  Maharudrayya  et  al.  [1]  numerically  ana¬ 
lyzed  the  effects  of  Reynolds  number,  aspect  ratio  and  curvature 
ratio  of  flow  channels  on  pressure  drop,  and  found  that  the 
overall  pressure  drop  reduces  with  the  number  of  the  serpen¬ 
tine  channels.  Also  by  numerical  analysis  of  fuel  cells  with 
serpentine  flow  channels,  Nguyen  et  al.  [2]  disclosed  that  the 
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activation  overpotential  mainly  occurs  under  the  ribs,  which 
implied  that  the  Ohmic  loss  under  the  ribs  is  lower  than  that 
under  the  channels.  In  addition,  it  was  indicated  that  the  temper¬ 
ature  at  the  interface  between  the  channels  and  diffuser  layer  is 
1-7  °C  higher  than  the  operating  temperature.  Using  a  commer¬ 
cial  codes,  FLUENT,  Dutta  et  al.  [3]  studied  three-dimensional 
mass  transport  process,  fuel  concentration  variation  and  their 
effects  on  electrochemical  reaction  in  PEM  fuel  cells  with  ser¬ 
pentine  flow  channels. 

Though  serpentine  flow  channel  design  promotes  fuel  cell 
performance,  it  has  penalty  of  higher  pressure  loss  as  a  result 
of  too  many  bends.  To  solve  this  problem,  Nguyen  [4]  pro¬ 
posed  interdigitated  flow  channel,  in  which  baffle  plates/blocks 
are  added  to  conventional  channels,  and  then  the  channels  are 
divided  into  two  parts,  inlet  and  outlet  channels.  Reactant  gases 
from  the  inlet  channels  are  forced  into  the  outlet  channels 
through  the  gas  diffusion  layer  (GDL).  With  the  assistance  of 
the  forced  convection  besides  diffusion,  gases  are  compelled  to 
pass  the  domain  under  the  ribs  in  a  shortest  way.  Um  and  Wang 
[5,6]  dealt  with  the  effects  of  oxygen  and  vapor  distributions  and 
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Nomenclature 

Ajgf. 

exchange  current  density  at  anode  side 

AjS;fc 

exchange  current  density  at  cathode  side 

a 

chemical  activity  of  water  vapor 

C 

mass  fraction 

Cs 

saturated  concentration  of  liquid  water 

CF 

quadratic  drag  factor 

c0 

connected  coefficient  between  source  term  and 
(00  -  05) 

D 

mass  diffusivity  (m2  s-1) 

Dh2 

hydrogen  diffusivity 

Do2 

oxygen  diffusivity 

D\\+ 

proton  diffusivity  in  membrane 

Dk,& ff 

effective  mass  diffusivity  for  the  kth  species 
(m2  s-1) 

^porous 

equivalent  surface  diameter  of  porous  media 

F 

Faraday  constant  (96,487  C  mol-1) 

H 

height  of  flow  field  (m) 

H+ 

hydrogen  ion 

H20 

water 

h 

local  current  density  in  the  y-direction  (Am-2) 

i 

current  density  (A  m-2) 

k 

current  density  at  anode  side  (Am-2) 

jc 

current  density  at  cathode  side  (Am-2) 

jo 

reference  current  density 

kc 

coefficient  of  water  vapor  condensation  rate  (s-1) 

h 

coefficient  of  water  vapor  evaporation  rate 

kp 

permeability  (m2) 

L 

length  of  flow  field  (m) 

M 

mole 

N 

number  of  baffles 

P 

pressure  (atm) 

P  sat 

saturated  water  vapor  pressure  (atm) 

P' 

perturbed  variation  of  pressure  in  control  volume 

Power 

power  (W  m-2) 

a 

source  term  of  variable  (0) 

R 

universal  gas  constant  (8.314  Jmol-1  K-1) 

s’ 

surface  area  (m2) 

Sc 

corrected  source  term  of  concentration  equation 

Si 

source  term  of  phase  potential  equation 

Sl 

source  term  of  liquid  water  effect 

Su 

source  term  of  momentum  equation  in  the 
v-direction 

sv 

source  term  of  momentum  equation  in  the 
y-direction 

Sw 

source  term  of  momentum  equation  in  the 
z-direction 

S 

volume  ratio  occupied  by  liquid  water 

T 

temperature  (K) 

t 

time  (s) 

u 

velocity  in  the  v-direction  (ms-1) 

V 

volume  (m3) 

Ycell 

operating  cell  pressure  (V) 

V 

velocity  in  the  y-direction  (ms-1) 

W 

width  of  flow  field  (m) 

w 

velocity  in  the  z-direction  (ms-1) 

X 

mole  fraction 

X 

velocity  in  the  v-direction  (ms-1) 

y 

velocity  in  the  y-direction  (ms-1) 

Zf 

spices  valence 

Z 

velocity  in  the  z-direction  (ms-1) 

Greek  letters 

£ 

porosity 

£eff 

effective  porosity 

A 

water  content  in  membrane 

V 

kinematic  viscosity  (m2  s-1) 

P 

density  (kgm-3) 

0 

over-potential  (V) 

electric  conductivity  in  membrane 

0 

dependent  variables 

r<p 

exchange  coefficient 

0 

phase  potential  function 

Ax 

height  tapered  ratio 

X 

width  tapered  ratio 

Superscripts 

ref 

reference  value 

T 

turtuosity  of  the  pore  in  the  porous  medium 

Subscripts 

a 

anode 

CL 

catalyst  layer 

c 

cathode 

channel 

channel 

eff 

effective 

GDL 

gas  diffuser  layer 

H+ 

hydrogen  ion 

h2 

hydrogen 

h2o 

water 

inlet 

entrance 

k 

the  kth  species  of  the  mixture 

Mem 

membrane 

02 

oxygen 

outlet 

outlet 

porous 

porous  medium 

sat 

saturation 

total 

total 

in  the  x- direction 

A 

in  the  y-direction 

z 

in  the  z-direction 

fuel  transport  on  electrochemical  reaction  in  three  dimensional 
fuel  cells  with  straight  channels  and  interdigitated  channels. 

Kazim  et  al.  [7]  and  Yi  and  Nguyen  [8]  analyzed  the  effects  of 
various  flow  channel  designs  on  PEM  fuel  cell  performance,  and 
concluded  that  interdigitated  design  is  superior  in  accelerating 
chemical  reaction  and  improving  cell  performance.  Hu  et  al. 
[9,10]  considered  in  detail  the  effects  of  interdigitated  design  on 
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cell  performance.  Their  results  showed  that,  relative  to  straight 
flow  channel,  interdigitated  design  in  the  cathode  side  causes 
higher  oxygen  concentration  and  lower  liquid  water  content  in 
the  diffuser  layer,  but  its  cell  performance  is  worse  if  the  cathode 
fuel  gas  supplied  without  humidification. 

Two  main  reasons  cause  water  molecules  appearing  within 
the  membrane.  One  is  the  electro-osmotic  drag  by  potential  gra¬ 
dient  and  the  other  is  the  back  diffusion  by  water  concentration 
gradient  across  the  cell.  When  these  two  kinds  of  movement  have 
the  same  rate,  the  water  in  the  membrane  lies  at  the  equilibrium 
state.  Okada  et  al.  [11]  developed  a  simple  model  to  investigate 
water  transport  in  the  membrane,  and  their  computational  results 
indicated  that  the  water  concentration  profile  is  influenced  by  the 
thickness  and  humidity  of  the  membrane  and  the  current  density 
of  the  cell.  Bussel  et  al.  [12]  divided  a  MEA  into  tens  of  units 
to  establish  a  two-dimensional  dynamical  model,  and  studied 
the  water  formation,  electro-osmosis  and  diffusion  process  in 
the  PEM,  and  their  effects  on  cell  performance  under  various 
operating  conditions.  Mazumder  and  Cole  [13,14]  demonstrated 
the  differences  between  two-dimensional  and  three-dimensional 
numerical  results,  and  employed  a  simplified  model  to  discuss 
the  effects  of  liquid  water  on  cell  performance.  Recently,  Yan 
and  co-workers  presented  a  series  of  experimental  and  numeri¬ 
cal  studies  about  the  flow  field  designs  on  the  cell  performance 
of  PEM  fuel  cells  [15-19]. 

As  mentioned  above,  up  to  now  the  flow  channel  designs 
of  PEM  fuel  cells  appeared  in  the  open  literature  are  limited  to 
straight,  serpentine,  Z-type  and  interdigitated  designs.  Basically, 
straight  channel  designs  have  lower  capability  in  cell  perfor¬ 
mance  and  liquid  water  removal,  though  their  pressure  drop  is 
lower.  Serpentine  designs  have  too  large  pressure  drop  due  to 
channel  bends.  Interdigitated  designs  also  have  very  high  pres¬ 
sure  drop  in  the  presence  of  channel  blockage.  Therefore,  based 
on  straight  flow  channels,  a  new  style  of  flow  field,  i.e.,  tapered 
flow  channel,  is  proposed  to  enhance  fuel  utilization  efficiency  at 
high  current  density.  In  this  work,  straight  flow  channels  of  taper 
ratio  in  height  or  width  are  proposed  in  order  to  improve  the  effi¬ 
ciency  of  fuel  transport,  fuel  utilization,  and  cell  performance 
in  PEM  fuel  cells.  For  cells  with  this  class  of  tapered  channel 
designs,  the  effects  of  flow  area  variation  along  the  flow  channel 
on  the  fuel  flow  velocity,  the  fuel  distribution,  the  fuel  transport 
and  utilization,  the  capability  of  the  liquid  water  removal,  and 
the  cell  performance  are  explored  numerically. 

2.  Analysis 


Table  1 

The  parametric  values  used  in  this  work 

Length  of  flow  channel  (L)  (mm) 

141 

Height  of  channel  (//channel)  (mm) 

1 

Width  of  channel  (W)  (mm) 

1 

Thickness  of  GDL  (//gdl)  (mm) 

0.3 

Thickness  of  CL  (//cl)  (mm) 

0.01 

Thickness  of  membrane  (//Mem)  (mm) 

0.035 

Porosity  of  GDL  (£gdl) 

0.4 

Porosity  of  CL  (scl) 

0.4 

Porosity  of  membrane  (£Mem) 

0.28 

Permeability  of  GDL  (£p,gdl)  (m2) 

1.76  x  KT11 

Permeability  of  CL  (£p,cl)  (m2) 

1.76  x  KT11 

Permeability  of  membran  (kp,Mem)  (m2) 

1.8  x  10“18 

Electric  conductivity  of  GDL  (£2-1  m-1) 

53 

Electric  conductivity  of  CL  (£2-1  m-1) 

53 

Reference  current  density  in  anode  side  (Ajj^)  (Am-3) 

9.23  x  108 

Reference  current  density  in  cathode  side  (AjJ}ctc)  (Am-3) 

1.05  x  106 

(3)  the  system  is  isothermal; 

(4)  the  flow  is  laminar; 

(5)  the  fluid  is  incompressible; 

(6)  the  thermal  properties  are  constant;  and 

(7)  the  porous  layers  such  as  GDL,  catalyst  layer  and  PEM  are 
isotropic. 


To  reduce  the  time  cost  of  three-dimensional  cell  computa¬ 
tion,  small  scale  fuel  cell  models  were  adopted  by  most  asso¬ 
ciated  literature.  However,  due  to  their  short  flow  channels,  the 
results  hardly  agreed  with  the  general  performance  of  the  cells 
at  a  large  scale.  This  work  aims  to  establish  a  relatively  large 
fuel  cell  model  with  the  area  of  14.1  x  14.1  cm1 2.  The  geometri¬ 
cal  dimensions  and  physical  properties  of  its  flow  channels  are 
listed  in  Table  1 .  Considering  the  symmetry  nature  of  the  fuel 
cell  configuration  as  shown  in  Fig.  1,  only  the  essential  part  of 
the  fuel  cell  model  is  taken  as  the  computational  domain  for 
saving  the  computational  time.  The  schematic  diagram  of  the 
three-dimensional  tapered  flow  channel  designs  is  depicted  in 
Fig.  2.  Two  parameters,  height  taper  ratio  Ax  and  width  taper 
ratio  Az 5  are  defined  as  follows: 


//outlet  . 

- ,  Az 

u  ’  4 

//inlet 


^outlet 

Wmlet 


where  //outlet,  //inlet,  ^outlet,  and  Winiet  are  the  outlet  height,  inlet 
height,  outlet  width,  and  inlet  width,  respectively. 

The  governing  equations  under  the  assumptions  mentioned 
above  can  be  depicted  as  follows: 


This  work  establishes  a  three-dimensional  PEM  fuel  cell 
model  to  analyze  its  electrochemical  reaction  and  transport  phe¬ 
nomena  of  reactants  in  anode  catalyst  layer,  PEM,  and  cathode 
catalyst  layer  in  the  MEA,  where  the  reaction  occurs.  In  other 
regions  such  as  the  GDL,  no  chemical  reaction  takes  place. 
Therefore,  the  governing  equations  of  this  model  include  mass, 
momentum,  species  and  electrical  potential  conversation  equa¬ 
tions.  For  simplicity,  the  following  assumptions  are  made: 

(1)  the  cell  system  is  three-dimensional  and  steady; 

(2)  the  inlet  fuel  is  ideal  gas; 


continuity  equation: 


du  dv  dw 
—  T  —  T  —  =  0. 
dx  dy  d  z 

momentum  equation  (v): 

du  du  du 
£eff  (  u  —  +  v—  +  w  — 

dz 


dx  '  dy 

£eff  3  P 


f)  dx 


d2u  d  2u  d2u 

+  U£eff| 


+  5 


u 


(3) 
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Fig.  1.  Schematic  diagram  of  the  three  dimensional  model. 


momentum  equation  (y): 


£eff 


dv  dv  dv 
u  - — b  v - b  w 


dx  dy 
£eff  dP 


3  z 


P  3  y 


d2v  d2v  d2v 

+  V£eff  (  + 


dx 2  dy 2  3  z2 


momentum  equation  (z): 


£eff 


3  w  dw  dw 

Li  - - b  V - b  w 


dx  '  dy 
£eff  3  P 


p  3  z 


+  D£eff 


3  z 
d2W 


+  s 


V 


d2w  3  2w 

dx2  +  dy2  +  dz 2 


w 


(4) 


(5) 


In  various  layers,  Su,  Sv  and  Sw  denote  different  source 
terms,  whose  details  are  listed  in  Table  2.  In  PEM,  Su ,  Sv  and 
Sw  embody  the  proton  transport  and  the  liquid  water  dragged 
by  electro-osmosis.  In  Table  1,  £eff  is  the  effective  porosity  of 
porous  material,  Cp  the  quadratic  drag  factor,  and  Zf  the  ion 
valence  in  PEM.  In  addition,  ff  is  the  effective  diffusion 
coefficient  of  the  Mi  speices,  which  is  calculated  by  the  B rug- 
german  equation  [20],  and  kv  is  the  permeability  in  porous 
material,  and  its  dependence  of  porosity  can  be  described  by 


Blake-Kozeny  equation  [21], 


/ d 2  \ 

1  ^porous  \ 

i 

00 

oo 

_ 1 

V  150  ) 

1 

(N 

CO 

1 

H 

_ 1 

where  dp0rous  is  the  equivalent  area  of  porous  material  and 
expressed  in  the  from  of 


d 


porous 


6Y 

w  y  porous 
S' 

^porous 


Species  equation  (Mi  species): 


where  Q  is  the  concentration  of  Mi  species,  and  Sc  is  the  source 
term  in  the  catalyst  layers  with  the  correction  by  chemical  reac¬ 
tion.  It  is  of  the  form  -jJ2FCtotSiU  for  hydrogen,  -7c/4FCtotai,c 
for  oxygen,  and  —jc/2FC{ Gtal,c  for  vapor.  In  these  expressions, 
7a  and  jc ,  respectively,  denote  the  current  densities  on  anode 
side  and  cathode  side,  and  can  be  calculated  by  the  following 


Table  2 

Source  terms  appeared  in  the  governing  equations 
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equations: 


k  =  AjJ) 


ref 


<%) 


Q(aaF/RT)r] 


1 

Q(acF/RT)r) 


Jc  =  Ajg 


ref 
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1 
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Besides  the  continuity,  momentum,  and  energy  equations, 
the  phase  potential  equation  for  potential  and  current  profile  is 


(10) 


where  &  is  the  phase  potential  function,  and  crm  is  the  mem¬ 
brane  conductivity  and  is  expressed  as  below  according  to 
Springer  et  al.  [22]  and  Nguyen  and  White  [23]. 


crm(T)  =  exp 


(11) 


where  crj^f  is  the  reference  conductivity  of  membrane,  and  is 
expressed  as  follows: 

o2f  =  0.005 139A  -  0.00326  (12) 


0.043  +  17.81a  —  39.85a2  +  36.0a3  0  <  a  <  1 

“  “  (13) 

14  +  1.4(a  —  1)  1  <  a  <  3 


in  which  a  denotes  the  vapor  activity  and  A  the  water  content 
in  PEM,  Sj  in  Eq.  (10)  is  the  electricity  source,  which  is  zero  in 
the  PEM  without  electrochemical  reaction  and  is  — ya  or  —jc  on 
the  cathode  or  anode  side.  The  relationship  between  the  phase 
potential  0  and  the  current  density  i  is  expressed  as  below: 


lx  —  OTi 


m 


ly  ~ 


—a, 


m 


d& 
dx 
d& 
3 y 
90 


(14) 

(15) 


E  —  Chn 


dy 


(16) 


Thereby,  the  phase  potential  equation  can  be  reduced  to  the  fol¬ 
lowing  equation: 


dix  diy  diz 

—  +  —  +  —  =  5,  (17) 

dx  dy  9z  J 

As  another  important  issue  of  this  model,  the  effects  of  liq¬ 
uid  water  are  simplified.  Although  there  is  no  two-phase  flow 
behavior  included  in  the  model,  the  blockage  effects  with  the 
presence  of  liquid  water  in  the  pores  of  porous  layers  such  as 
the  GDL,  catalyst  layer  and  PEM,  are  taken  into  account.  The 
source  term  5l  in  the  species  equations  is  expressed  to  be  the 
water  amount  of  condensation  or  evaporation,  and  is  evaluated 
as: 


Mn2okc  efpj^°  (^H2°  “  ^sat),  ifftfcO  >  ^sat 
ke£effS(PSat  -  ^H2o),  if  ^H20  <  ^sat 

(18) 
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Cathode  bipolar  layer 
Cathode  diffuser  layer 
Cathode  catalyst  layer 
Membrane 
Anode  catalyst  layer 
Anode  diffuser  layer 
Anode  bipolar  plate 


Fig.  2.  Schematic  diagrams  of  the  tapered  channel  designs:  (a)  flow  channel  tapered  in  height;  (b)  flow  channel  tapered  in  width. 


where  M  is  the  molar  weight,  kc  and  kQ  the  vapor  condensation 
and  evaporation  constant,  and  Psa t  the  vapor  saturated  pressure, 
which  is  calculated  as  follows  [16-18]: 

p  =  i0_2-1794+0-02953r_9-1837xl0_5r2+L4454xl0_7r3  (19) 


where  s  denotes  the  volume  ratio  of  liquid  water  to  pore  of  the 
porous  material,  and  that  is, 


volume  of  liquid  water 
volume  of  pore 


(20) 


Liquid  water  affects  the  effective  porosity  of  porous  layers, 
which  is  expressed  as  below: 


£eff  =  e(l  -  s) 


(21) 


2.1.  Boundary  conditions 


The  boundary  conditions  of  the  present  computational  model 
contain  those  at  inlet,  outlet,  wall,  and  interfaces  between  differ¬ 
ent  layers.  For  comparison  with  the  experiments,  the  inlet  flow 
rates  in  cathode  and  anode  flow  fields  are,  respectively,  specified 
to  be  2000  and  4000  cm3  min-1 ,  and  the  other  specifications  are 
listed  in  Table  3.  At  the  outlet,  fully  developed  flow  is  assumed 
and  the  boundary  conditions  of  velocity  and  concentration  can 
be  expressed  as: 


dv 

dy 


=  w  = 


(22) 


At  the  interface  between  the  flow  channels  and  bipolar  plates, 
the  velocity  and  mass  flux  are  zero. 


dck 

u  =  v  =  w  =  - =  0 

3  z 


(23) 


At  the  interface  between  the  flow  fields  and  the  GDLs,  the 
velocity,  mass  flux,  species  concentration  and  species  flux  are 
all  continuous. 


£eff,GDL 


dv 

dx 


*=^GDL 


dv 

dx 


£eff,GDL 


X — ^channel 


dw 

dx 


X=XqGL 


dw 

dx 


Vx=XGdl  ~  vx=X 


channel 


(24) 


X — ^channel 


fceff,GDL 


dx 


x=Xqdl 


dx 


k>X — ^GDL  ^ k,X — ^channel 


(25) 


Table  3 

The  inlet  conditions  of  fuel  gases  used  in  this  work 


Anode  side 

Cathode  side 

Inlet  velocity  (ms-1) 

0.475 

0.95 

Fuel 

h2,  h2o 

o2,  h2o,  n2 

Relative  humidity  (%) 

100 

100 

Temperature  (°C) 

50 

50 

U  = 


W.-M.  Yan  et  al.  /  Journal  of  Power  Sources  161  (2006)  907-919 


913 


At  the  interface  between  the  GDLs  and  the  catalyst  layers, 
the  velocity  and  concentration  conditions  are  similar  to  those  at 
the  interface  between  the  flow  fields  and  the  GDLs. 


£eff,CL 


£eff,CL 


dv 

dx 

dw 

dx 


=  £eff,GDL 


X=XCL 


dv 

dx 

dw 


x=Xqdl 


=  £eff,GDL 


x=XqL 


dx 


*=^GDL 


Vx=xCL  Vx=Xqdl  ’  wx=Xcl  wx=Xgdl 


(26) 


dCk 

£eff,CL  “ 
dx 


x=Xql 
Ck,x=XcL  =  Ck,X=XQY>L 


=  ^eff,GDL 


dx 


X=Xqd1 


(27) 


Also,  at  the  interfaces  between  the  catalyst  layers  and  PEM, 
as  gases  cannot  permeate  PEM,  their  velocity,  mass  flux,  species 
concentration  and  species  flux  are  all  zero. 


where  t  is  the  time,  0  general  dependent  variable,  such  as  veloc¬ 
ity  components,  concentration,  exchange  coefficient  such  as 
viscosity,  diffusion  coefficient,  source  term  such  as  pressure 
gradient,  u  velocity  vector  and  p  is  the  density. 

When  Eq.  (32)  is  integrated  over  a  control  volume,  a  finite- 
volume  equation  can  be  resulted  in  as 

<2p0p  =  <3E0E  +  tfw0W  +  ^N0N  +  «S0S  +  ^H0H  +  ^L0L  +  b 

(33) 

where  0p  is  the  value  of  variable  0  at  node  P  of  control  volume, 
0E,  . . . ,  0l  value  of  variable  at  node  of  neighboring  control  vol¬ 
ume,  <2p,  . . . ,  <2l  coefficients  in  discretized  equations  and  b  is 
the  expression  of  source  for  the  discretized  equation  of  variable 
4>. 

The  Semi  Implicit  Method  for  Pressure-Linked  Equation 
(SIMPLE)  algorithm,  developed  by  Patankar  [24],  is  employed 
to  solve  the  governing  equations. 

3.  Results  and  discussion 


dv 

£eff,CL  — 

dX  *=*CL 

Vx=XCL  =  0 


=  0,  £eff,CL 


dw 

dx 


=  0, 


*=^CL 


(28) 


S  eft,  CL 


dCk 

dx 


=  o,  Ck,x=x  CL  =  0 


x=Xql 


(29) 


Since  liquid  water  can  transport  through  the  PEM  as  a  result 
of  the  influences  of  proton  transport  and  electro-osmosis,  its 
velocity,  mass  flux,  liquid  water  saturation  degree  Cs  and  its 
flux  is  continuous  at  the  interfaces  between  the  catalyst  layers 
and  the  PEM. 


£eff,Mem 


£eff,Mem 


dv 

dx 

dw 

dx 


=  £eff,CL 


x=Xm 


em 


=  £eff,CL 


X — Xvicm 

Vx=xMcm  =  Vx=xCh, 


dv 

dx 

dw 

dx 


x=Xcl 


x=Xql 
WX=X  Mem  =  WX=Xql 


(30) 


The  validation  of  the  numerical  results  is  performed  by 
comparing  the  present  predictions  with  previous  experimental 
results.  Fig.  3  shows  the  comparison  of  the  numerical  and  experi¬ 
mental  results  of  the  cell  performance,  i.e.,  the  operation  voltage 
Vceii  versus  the  current  density  of  fuel  cells  I  (A  m-2),  at  the 
operating  conditions  of  fully  humidified  fuel  gases  with  flow 
rate  of  4000  cm3  min-1 ,  pressure  of  1  atm,  and  humidifying  and 
operating  temperatures  of  50  °C.  The  solid  circles  denote  the 
experimental  results  of  Cheng  [25]  and  the  curve  is  the  present 
numerical  prediction.  The  numerical  results  show  a  few  devi¬ 
ations  from  the  measurements  at  low  and  high  current  density 
regimes.  Since  in  the  numerical  model  the  presence  of  the  liquid 
water  is  assumed  to  cause  only  the  effects  of  blocking  pores  and 
hindering  fuel  transport.  This  is  different  from  the  complicated 
two-phase  transport  phenomena  in  a  real  case.  At  high  current, 
the  numerical  analysis  with  a  simplified  model  underestimates 
the  liquid  water  effects  and  the  model  is  not  appropriate  to  simu¬ 
late  the  real  dead-zone  condition.  Nevertheless,  the  discrepancy 
of  the  cell  performance  between  the  numerical  predictions  and 


£cff,  Mem 


dCs 

dx 


=  £eff,CL 


X — Vfvlcm 

Cs7x=V|vlem  =  Cs,X=XCL 


dCs 

dx 


x=Xcl 


(31) 


The  boundary  conditions  of  potential  at  the  interface  between 
the  anode  catalyst  layers  and  the  PEM  is  0  =  0,  and  that  at  the 
interface  between  the  cathode  catalyst  layer  and  the  PEM  can 
be  expressed  as  d@/dx  =  d@/d z  =  0. 

The  general  convection-diffusion  equation  can  be  expressed 
in  a  conservative  form  as  follows: 


— (f>4>)  +  V  •  (/ ou<p 
ot 


E<pV4>)  =  S<p 


(32) 


Fig.  3.  Comparison  of  the  numerical  predictions  with  the  experimental  results. 
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Fig.  4.  Effects  of  tapered  channel  design  in  2-D  and  3-D  models  on  the  I-V 
curves. 

the  experimental  results  lies  within  a  reasonable  range.  Thereby 
it  can  be  deduced  that  the  present  numerical  method  is  appro¬ 
priately  valid  in  simulation  of  the  three-dimensional  fuel  cells 
with  the  presence  of  liquid  water. 

Using  this  numerical  model,  the  tapered  flow  channel  design 
is  analyzed  with  the  emphasis  on  the  effects  of  the  height  taper 
ratio  and  width  taper  ratio  on  cell  performance.  Fig.  4  presents 
the  effects  of  height  taper  ratio  on  the  two-dimensional  (2-D) 
and  three-dimensional  (3-D)  predictions  of  the  cell  performance. 


Fig.  5.  Effects  of  height  taper  ratio  Ax  on  the  cell  performance,  (a)  I-V  curves; 
(b)  I-P  curves. 


Fig.  6.  Effects  of  width  taper  ratio  Az  on  the  cell  performance,  (a)  I-V  curves; 
(b)  I-P  curves. 

The  solid  and  dashed  lines,  respectively,  denote  the  results  of  3- 
D  and  2-D  models.  It  is  disclosed  in  Fig.  4  that,  for  the  voltage 
less  than  0.6  V,  the  current  density  predicted  using  the  2-D  model 
is  about  two  times  of  that  using  the  3-D  model,  because  the  com¬ 
putational  domain  of  the  2-D  model  is  equivalently  the  central 
section  of  the  3-D  model  or  section  A  in  Fig.  1,  where  the  cur¬ 
rent  density  is  the  maximum  in  the  3-D  model.  The  diffusion 
path  of  the  fuel  under  the  ribs  is  longer  than  that  under  the  flow 
channels,  and  part  of  fuel  under  the  flow  channels  has  partici¬ 
pated  in  the  reaction,  so  the  concentration  of  the  fuel  under  the 
ribs  is  lower.  Since  the  concentration  of  the  fuel  under  the  ribs 
is  low  and  the  liquid  water  is  difficult  to  be  removed,  the  current 
density  lowers  down,  and  thus  the  current  density  calculated  by 
the  3-D  model  are  less  than  that  by  the  2-D  model. 

Figs.  5  and  6  present  the  effects  of  the  height  taper  ratio  Ax 
and  the  width  taper  ratio  Az  on  cell  performance,  and  Fig.  5 (a 
and  b)  denotes  the  results  of  I-V ceii  (current  density-voltage) 
and  I-P  (current  density-power),  respectively.  At  high  voltage 
(Vceii  >  0.7  V),  the  tapered  channel  design  has  little  effects  on 
current  density  for  low  electrochemical  reaction  rate.  As  the 
voltage  decreases  and  the  current  density  increases,  the  fuel 
amount  needed  rises  up  for  more  active  electrochemical  reac¬ 
tion,  and  the  limiting  condition  due  to  lack  of  fuel  occurs  as  the 
voltage  lies  at  about  0.5  V.  For  the  lower  height  taper  ratios  in 
Fig.  5(a),  fuel  transport  into  GDL  can  be  enhanced  in  the  region 
near  the  end  or  the  downstream  portion  of  the  flow  channels,  so 
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Fig.  7.  Effects  of  height  taper  ratio  Ax  on  the  local  current  density  distributions 
at  various  cross  sections  of  the  membrane,  (a)  V^n  =  0.7  V  and  (b)  Vceii  =  0.3  V. 


Y(m) 


Fig.  8.  Effects  of  width  taper  ratio  Az  on  the  local  current  density  distributions 
at  various  cross  sections  of  the  membrane,  (a)  Vcen  =  0.7  V  and  (b)  Vceii  =  0.3  V. 


the  fuel  amount  increases  and  the  occurrence  of  the  fuel  mass 
transport  limit  can  be  delayed.  In  addition,  at  a  lower  current 
density  (/<  4000  A  m-2),  the  height  taper  ratio  has  little  effects 
on  the  output  power  for  slow  electrochemical  reaction,  while 
at  higher  current  density  (/>  4000  A  m~2),  a  lower  height  taper 
ratio  results  in  a  higher  output  power.  With  increasing  width 
taper  ratio,  the  contact  area  between  fuel  and  GDL  rises  for  the 
area  increment  of  the  flow  channels,  then  the  content  of  fuel  in 
the  GDL  becomes  higher  and,  subsequently,  the  current  den¬ 
sity  also  rises.  Fig.  6(b)  shows  the  growth  of  the  output  power 
with  the  increasing  width  taper  ratio  corresponding  to  the  cases 
presented  in  Fig.  6(a). 

Figs.  7  and  8  present  the  effects  of  the  height  and  width  taper 
ratios  on  local  current  distributions  at  various  sections.  Along 
the  fuel  channels,  the  fuel  amount  drops  due  to  consumption 
in  reaction,  while  the  liquid  water  as  a  product  of  the  reaction 
accumulates.  As  a  result,  at  the  end  of  the  flow  channels,  the 
fuel  transport  rate  lowers  down,  and  the  local  current  density 
decreases  along  the  channel  (F)  direction.  The  data  in  Fig.  7(a) 
show  that  the  local  current  density  increases  with  the  height 
taper  ratio  in  the  inlet  region  (Y<  0.02  m),  because  the  liquid 
water  accumulates  in  this  region  and  increases  the  resistance  of 
fuel  diffusion.  Along  the  channel  length,  more  fuel  in  the  flow 
channels  diffuses  to  the  GDL  and  catalyst  layer.  This  enhance¬ 
ment  is  caused  by  the  contraction  of  the  flow  area  of  the  fuel 
channels  with  its  height  tapered  along  the  channel  length.  As  a 


result,  the  lower  Ax  generates  better  current  density  distribution 
in  the  region  of  Y >  0.02  m.  In  the  inlet  region  of  the  flow  chan¬ 
nels  of  Ax  =  0.1,  the  fuel  diffusion  to  the  GDL  is  enhanced  due 
to  the  channel  taper  effect,  which  not  only  increase  the  current 
density  and  fuel  flux  but  also  benefit  the  liquid  water  removal.  At 
the  voltage  of  0.3  V  in  Fig.  7(b),  the  current  density  distribution 
is  similar  to  that  at  0.7  V.  However,  since  the  electrochemical 
reaction  at  0.3  V  is  very  drastic,  a  large  amount  of  fuel  reacts  in 
the  region  beneath  the  flow  channels.  It  leads  to  very  low  current 
density  under  the  ribs,  so  the  effects  of  height  taper  ratio  are  not 
evident  under  the  ribs. 

In  Fig.  8(a),  the  effects  of  the  width  taper  ratio  at  the  operat¬ 
ing  voltage  of  0.3  V  are  quite  obvious.  At  Az  >  1,  for  its  larger 
area  of  the  flow  channels,  the  fuel  velocity  is  lower  than  that  at 
Az  <  1 .  With  an  increase  in  the  channel  width,  the  fuel  has  more 
time  to  diffuse,  and  its  diffusion  path  under  the  ribs  decreases, 
so  the  difference  between  the  current  density  distributions  in 
the  regions  under  the  ribs  and  the  channels  becomes  smaller.  At 
the  largest  taper  ratio  Az=  1.8,  the  local  current  density  at  the 
outlet  is  almost  equal  to  that  under  the  flow  channels.  At  the 
operating  voltage  of  0.3  V,  the  results  shown  in  Fig.  8(b)  reveal 
that  the  local  current  density  drops  very  quickly  along  the  chan¬ 
nel  because  most  of  the  fuel  takes  part  in  the  electrochemical 
reaction.  For  low  width  taper  ratio,  the  fuel  and  the  associated 
reaction  concentrate  in  the  central  region.  As  a  result,  with  a 
reduction  in  width  taper  ratio,  the  local  current  density  under 
the  flow  channels  rises,  while  under  the  ribs  it  decreases  for 
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Fig.  9.  Effects  of  tapered  channel  design  on  the  velocity  distributions  along  the 
cathode  flow  channels:  (a)  tapered  in  height;  (b)  tapered  in  width. 

the  effects  of  diffusion  path.  In  the  case  of  the  flow  channels 
with  smaller  area,  the  local  current  density  under  the  flow  chan¬ 
nels  is  generally  higher,  while  under  the  ribs  the  current  density 
is  relatively  lower  for  the  concentrated  reaction  under  the  flow 
channels,  so  larger  width  taper  ratio  produces  higher  current 
density. 

The  fuel  velocity  in  the  flow  channels  affects  the  liquid  water 
removal  in  fuel  cells.  Fig.  9  presents  the  effects  of  tapered  chan¬ 
nel  on  the  fuel  velocity  at  the  center  of  the  channel.  The  effects 
of  height  and  width  tapered  channels  are  shown  in  Fig.  9 (a  and 
b),  in  which  the  continuous  and  dotted  lines  denote  the  results  at 
0.7  and  0.3  V,  respectively.  In  the  inlet  region  of  the  flow  chan¬ 
nels,  the  flow  develops  approaching  fully  developed  condition, 
and  the  center-line  velocity  usually  has  a  significant  increment. 
In  the  baseline  case  of  Ax=  1 ,  the  velocity  changes  little  and 
slightly  increases  along  the  channel  for  the  consumption  of  the 
fuel  and  the  decrease  of  the  fuel  density  on  the  cathode  side. 
For  the  design  of  height  tapered  flow  channels,  as  shown  in 
Fig.  9(a),  the  reduction  in  cross  section  area  accelerates  the  fuel 
flow  along  the  channels  like  the  flow  in  a  nozzle.  This  channel 
tapered  effect  becomes  more  remarkable  as  Ax  decreases.  In  the 
width  tapered  flow  channels  of  Az  <  1,  due  to  nozzle  effect  and 
the  fuel  consumption,  the  fuel  velocity  increases  along  the  chan¬ 
nel,  and  the  acceleration  is  more  evident  for  lower  Ax.  For  the 
flow  channels  of  Az  >  1 ,  the  center- line  velocity  drops  along  the 
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Fig.  10.  Effects  of  tapered  channel  design  on  the  oxygen  distributions  along  the 
cathode  flow  channels:  (a)  tapered  in  height;  (b)  tapered  in  width. 

flow  for  the  increasing  cross-sectional  area,  although  the  fuel 
consumption  has  somewhat  influences  in  increasing  the  veloc¬ 
ity.  In  the  case  of  a  lower  fuel  velocity,  the  liquid  water  removal 
caused  by  electrochemical  reaction  is  difficult,  but  the  fuel  can 
have  enough  time  for  diffusion  towards  the  catalyst  layer.  At  low 
voltage,  most  oxygen  participates  the  reaction,  which  results  in 
low  fuel  density.  Therefore,  the  local  velocity  is  higher  than  that 
at  a  higher  voltage.  Moreover,  this  situation  is  more  obvious  as 
Ax  decreases. 

Because  the  tapered  flow  channels  are  of  a  reduction  in  the 
channel  cross  sectional  area,  the  fuel  in  it  can  be  more  easily 
to  diffuse  to  the  porous  layers.  Fig.  10  shows  the  effects  of  the 
tapered  channels  on  fuel  distributions  (weight  concentration  of 
oxygen)  on  the  cathode.  Fig.  10(a  and  b)  displays  the  results 
of  the  height  and  width  tapered  flow  channels,  respectively.  At 
operating  potential  of  0.7  V,  for  height  tapered  channels,  the 
results  in  Fig.  10(a)  disclose  that  the  fuel  continuously  diffuses  to 
the  GDL  and  catalyst  layer  and,  then,  takes  part  in  the  reaction  in 
the  catalyst  layer  and  the  oxygen  concentration  decreases  along 
the  channel.  Whereas  Ax  has  only  little  effects  on  the  oxygen 
concentration  for  the  slow  electrochemical  reaction.  However, 
at  a  low  voltage,  0.3  V,  the  current  density  gets  higher  and  more 
oxygen  is  consumed.  Besides,  the  taper  effects  are  more  obvi¬ 
ous.  A  smaller  Ax  shows  more  evident  taper  effects  and  forces 
more  fuel  to  reach  catalyst  layer  and  react  there.  It  implies  that 
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Fig.  1 1 .  Effects  of  tapered  channel  design  on  the  liquid  water  distributions  along 
the  cathode  flow  channels:  (a)  tapered  in  height;  (b)  tapered  in  width. 

the  tapered  channels  increase  the  fuel  utilization  efficiency.  For 
flow  channels  with  large  Az,  as  shown  in  Fig.  10(b),  the  aver¬ 
age  oxygen  concentration  reduces  due  to  high  utilization  of  the 
fuel  oxygen.  It  is  also  observed  that  at  a  high  current  density, 
the  drastic  electrochemical  reaction  consumes  more  oxygen,  but 
some  oxygen  in  the  flow  channels  cannot  transport  to  the  catalyst 
layer.  Since  a  flow  channel  of  a  large  width  taper  ratio  have  rela¬ 
tively  larger  flow  area,  which  is  beneficial  to  the  fuel  utilization 
efficiency. 

Effects  of  the  tapered  channel  on  the  liquid  water  distribution 
in  the  cathode  are  shown  in  Fig.  1 1 ,  in  which  the  parameter  S  is  an 
index  denoting  the  degree  of  the  liquid  water  saturation.  The  liq¬ 
uid  water  content  increases  along  the  channel  as  a  consequence 
of  the  electrochemical  reaction.  Fig.  11  (a  and  b)  illustrates  the 
effects  of  height  taper  ratio  and  width  taper  ratio,  respectively. 
At  the  voltage  of  0.7  V  and  the  position  of  Y<  0.05  m,  the  liq¬ 
uid  water  distribution  is  similar  for  the  small  taper  effect.  In 
the  downstream  portion  of  the  channel,  as  shown  in  Fig.  9(a), 
the  velocity  rises  in  the  tapered  channel  and  more  liquid  water 
can  thus  be  removed.  As  a  result,  in  the  downstream  part  of  the 
channels  (Y >  0.05  m),  the  liquid  water  distribution  of  a  tapered 
channel  is  lower  than  that  of  the  regular  straight  channel  with 
Ax  =  1 .  As  Ax  =  0. 1  and  Y >  0.075  m,  the  liquid  water  removal  is 
performed  well  for  the  velocity  increasing  with  Y  as  a  result  of  the 
large  taper  effect.  It  is  shown  that  in  this  region  the  liquid  water 


content  rises  with  the  increase  in  Ax.  At  0.3  V,  the  more  drastic 
electrochemical  reaction  produces  more  liquid  water  than  that 
at  0.7  V,  and  the  region  with  evident  tapered  effect  expands  from 
Y >  0.05  m  to  Y >  0.025  m.  Near  the  outlet,  the  taper  effects  are 
stronger,  and  the  difference  of  liquid  water  distribution  is  more 
obvious.  In  Fig.  1 1(b),  at  0.7  V  and  in  the  region  of  Y<  0.05  m, 
the  width  of  the  flow  channels  changes  little,  and  the  liquid  water 
produced  in  the  catalyst  layer  has  not  permeate  back  to  the  flow 
channels,  so  the  width  taper  ratio  effects  on  liquid  water  distri¬ 
bution  are  not  apparent.  However,  along  the  channel,  the  liquid 
water  continuously  flows  in  the  channel  and  accumulates  in  the 
downstream  region.  For  a  tapered  channel  of  Az  <  1,  the  flow 
velocity  increases  along  the  channel,  and  the  liquid  water  is  not 
easy  to  accumulate  in  the  channel;  while  in  a  divergent  channel, 
Az>  1,  more  fuel  takes  part  in  reaction  and  more  liquid  water 
forms.  The  flow  velocity,  as  illustrated  in  Fig.  9(b),  decreases 
along  the  channel,  so  the  liquid  water  distribution  increases  with 
the  increasing  width  taper  ratio  of  the  flow  channel.  These  phe¬ 
nomena  are  more  obvious  at  0.3  V,  in  which  the  liquid  water 
distribution  significantly  increases  for  the  high  reaction  rate,  and 
it  rises  by  1-2  times  of  that  at  0.7  V.  Therefore,  a  smaller  width 
taper  ratio  leads  to  better  capability  of  liquid  water  removal. 

The  tapered  channel  design  is  beneficial  to  the  improvement 
in  fuel  diffusion  and  cell  performance,  and  it  also  alters  the  resis¬ 
tance  of  the  fuel  flow.  In  order  to  explore  the  pressure  loss  caused 
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Fig.  12.  Effects  of  tapered  channel  design  on  the  pressure  distributions  along 
the  cathode  flow  channels:  (a)  tapered  in  height;  (b)  tapered  in  width. 
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by  the  flow  area  contraction  in  tapered  channels,  the  local  pres¬ 
sure  losses  (the  difference  between  local  and  inlet  pressures)  at 
various  locations  of  the  center-line  of  the  cathode  flow  channels 
are  presented  in  Fig.  12.  As  illustrated  in  Fig.  12(a),  the  local 
pressure  loss  increases  with  channel  length.  At  Ax  =  0.5,  the  dif¬ 
ference  between  the  inlet  and  outlet  pressures  is  about  two  times 
of  that  at  Ax  =  1 .  At  Ax  =  0. 1,  the  pressure  loss  almost  rises  lin¬ 
early  for  Y <  0.08  m,  while  it  rises  in  a  nonlinear  manner  with  an 
increasing  taper  effect.  As  the  main  reason,  when  Ax  =  0.1,  the 
height  of  inlet  is  just  ten  percent  of  the  original  one,  and  both 
the  flow  resistance  and  pressure  loss  are  very  large.  As  shown 
in  Fig.  12(b),  the  design  with  smaller  width  tapered  ratio  of 
Az  <  1  can  increase  the  fuel  flow  velocity  and  improve  the  per¬ 
formance  of  liquid  water  removal,  but  the  resistance  of  fuel  flow 
also  increases;  the  design  with  larger  width  taper  ratio  Az>  1 
increases  the  diffusion  area  of  fuel,  but  there  accumulates  more 
liquid  water  in  the  flow  channels. 

4.  Conclusions 

In  the  present  paper,  the  relatively  novel  channel  design  of 
variable  cross-section  has  proposed  and  the  effects  of  the  taper 
ratio,  either  in  height  or  width,  have  been  analyzed  in  details. 
With  this  class  of  fuel  channels,  the  present  numerical  analysis 
has  demonstrated  that  the  flow  area  variation  along  the  flow 
channel  leads  to  noticeable  influences  on  the  fuel  velocity  in 
flow  channel,  the  capability  of  the  liquid  water  removal,  the  fuel 
transport  and  utilization,  and  the  cell  performance.  Based  on  the 
present  predictions  and  analysis,  the  conclusions  can  be  derived 
as  follows: 

1 .  For  the  cathode  channels  tapered  in  height,  the  reduction  in 
outlet  height  is  beneficial  to  the  fuel  utilization  efficiency  and 
the  improvement  in  the  cell  performance. 

2.  For  the  tapered  channel  with  the  width  of  flow  channel 
enlarged  along  the  flow,  the  contact  area  of  the  interface 
between  channels  and  GDL  increases,  which  in  turn,  results 
in  an  increase  in  the  transport  rate  of  fuel  in  the  outlet.  There¬ 
fore,  the  cell  performance  is  improved. 

3.  The  flow  channels  tapered  in  height  can  enhance  the  fuel 
velocity  along  the  flow  channel,  and  improve  the  liquid  water 
removal,  whose  rate  increases  with  the  decrease  in  height 
taper  ratio.  For  the  design  of  width  tapered  flow  channels,  the 
rise  of  the  width  decelerates  the  flow  and  weakens  the  capa¬ 
bility  of  liquid  water  removal,  while  it  increases  the  reaction 
efficiency  and  cell  performance. 

4.  At  the  height  taper  ratio  Ax  of  0.1,  the  design  with  a  large 
width  tapered  ratio  weakens  the  taper  effect.  Therefore,  the 
cell  performance,  the  fuel  velocity  and  the  efficiency  of  liq¬ 
uid  water  removal  are  all  reduced  as  the  width  taper  ratio 
increases. 

5.  If  the  power  loss  due  to  pressure  drop  is  not  considered,  the 
performance  of  the  fuel  cell  with  the  tapered  flow  channels 
is  consistently  improved  with  height  taper  ratios  decreased 
and  width  taper  ratio  increased.  With  the  pressure  loss  con¬ 
sidered,  however,  the  best  performance  can  be  obtained  at 
the  height  taper  ratio  ( Ax )  of  0.5  and  the  width  taper  ratio 


( Az )  of  1.8  among  the  taper  ratios  studied  in  the  present 
work. 
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